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INTRODUCTION 


The purpose of this project vs s to develop optimum techniques 
for training and consulting with users and potential users of 
Landsat digital computer technology. Specific goals were to: 

1) Develop technical schedules and consulting techniques for 
Integrating Landsat classifications with geo-based information 
systems (CIS) through cooperation with the U.S. Forest Service, 
McCloud Ranger District, data base integration project. 

2) Develop appropriate methodologies for expanding the information 
displayed by an unsupervised Landsat classification of the State 
of California for a limited geographical region. This investi- 
gation was to be done in cooperation with the California 
Department of Forestry (CDF) in a northern California county. 

3) Produce a Semi-annual and a Final Report documenting significant 
results and conclusions. 

The first goal, integration of a GIS, could not be accomplished since 
the Forest Service decided against incorporating Landsat data into 
their GIS at this time. It was therefore decided to alter the grant 
through mutual agreement between the Technical Officer and the 
Principal Investigator. The first goal was changed to: developing 
Landsat training techniques through the support of a VICAR Software 
Workshop at the NASA Ames Research Center. The material developed 
for this workshop and a short summary statement of results are pre- 
sented in Section One of this report. 

The second goal, Landsat information extraction, was accomplished as 
planned and a description of the findings is included as Section Two. 
The semi-annual report has been submitted and this document is the 
Fiual Report, fulfilling goal number three. 



SECTION ONE 


VICAR Workshop 


Project Overview 

The Principal Investigator developed and tested a teaching nodule 
on image classification procedures using the VICAR computer software 
package. The module was developed to optimize the training benefits 
for State users of the VICAR programs. The module was field tested 
at a NASA sponsored VICAR users workshop in the Fall of 1980. The 
module is organized into three basic sections: (1) Lecture on Image 

Classification using VICAR, (2) Exercise on classification with VICAR 
and (3) Discussion of results. 
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13. Unacceptable Conditions 
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IV. Computer Classification 

A. Gcrora? problem - Subdivide spectral space such that: 

1. AH spectral clauses ore separable . 

2. A.'l spectral classes leprcse; t ore irfonv-ica does. 

B. O', norr.3 Procedure 

1. Define or eh spectral clans ' *or c.-vory char.rel) by: 

a) Kern refit ctrr.cc value. 
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4. Evaluate results. 

5. Iterate if accuracy is poor. 
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C. Problems: 
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or trainin'; area. !«o way of knowing this. 

General con Ci.pt subdivides spectral space into information classes 


r.ot spectral classes. 

Ur.suporviced Classification or Clustering 

A . Approach 

1. Select number of spectral classes desired. 

2. Use cJucterinc routine to establish spectra?, classes, a fi:u\l 
stat file. 

J.. Determine the info nation class associated with each spectral 
class defined. 

3. When used: resource inforuatior. classes represented by severa 1 

spectral classes. 

C . Problems 

1. Lir.its on the number of pixels clustered at one time. 

2. Pixels containin':; two or c inform !: ion classes arc sample'’ 

•.-•inch r.:ay distor* statistic:*. 

J,. Limited control of v)*,: t el.. anon are created. 
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G lided Cluster ny 
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C. Pool a: nilar classes or delete clnsroe as necessary- to create 

r. tnt Tile. 

' . Kozinara livelihood classification. 

71. Aevnatacoe 

1. Control over tic* clustering; process as a linitct number o_ 

s. 'octroi c locoes are recreated within a training area, 

2. Pixels cor.tair.rns nony info mat ion classes are relocated to 
n separate class which can bo deleted . 

3. C3usterir.s finds spectral variation within inioroation c-... . o. . 
h. P* csible to concentrate on separating inforuation classes \;hic.: 

arc similar spectrally. 

C . Problems 

1. Spectral classes r.ot occv-risf, in trnininG fields are rot con- 
cidcrod, causing ciesclocsificntior.. 

2. Difficult to establish trsininp areas contnir.inG enough pixels 
(several spectral classes per training area). 
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be tailed Discussion of Work Flow 

K 
t 

1 • Select Study Area 

A. Site 

) 

The selection of an area for computer classification is influenced 
by several factors: 

1) Management unit or region of responsibility 

2) Natural boundaries of ecosystems, 6uch as watersheds 

3) For experimental work, the degree to which hind-use categories of 
interest are represented on the ground. 

The instructors have selected a study area for the exercises today. 
It is a coastal region in Northwestern Columbia, South America. The 
region was selected because it illustrates the international scope of 
the landsut system, contains several general categories of land-use, 
and we happened to have the Landsat CCT sitting on the shelf. 

B. Land Cover Information Categories 

The selection of information categories is the beginning step in 
classifying landsat data. Information category selection can be thought 
of us goal setting in the planning process. Dcciuiono must be irad<; on 
what land cover types are needed for the inventory and what landsat 
spectral signatures are apt to represent those cover types. 

Examining air photos of the study area will give an indication of 
the land covor categories present. We have some space photography from 
Shy lab and some high resolution KBV imagery of the study area foi you 
to use in this exercise. An enhanced londsnt, false color composite 
of the study area is also available. You should oelect approx irately 
3-b categories of land-use which you bolievc to be present in the study 
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JI. Select 6 Training Are an 

Training areas should be selected from grey scale computer print- 
out of the 6tudy area. These areas should be representative canplcs 
of the information categories you have selected and contain at leant 
100 pixels. Avoid border pixels. Coapare the grey scale print-out 
with the photos and false color composite to locate training areas 
for about 3 of your recource categories that appear very homogeneous. 
Outline, in pencil, training areas for 3 of your resource categories 
on the grey scale print-out, in rectangular fora. VICAR con accept 
irregular shapes through the "VERT" format but we will stay with rec- 
tangular today. Record the otarting line, starting soaplc, number of 
lines and number of samples .or each training area. We would like you 
to select 3 training areas. This means you have room for 2 renourcc 
categories to be represented by 2 training areas esch. The remaining 
resource categories are to be represented by one training area. In 
an actual application run, one would want several training areas to be 
specified for each reuourco category of interest. Thio allows fox* a 
more representative sample to be selected. 

Nov, you will probably find thut at least one of your resource 
categories did not fit nicely into a rectangular, homogeneous area. 

This is u very roiiscn problem in landaat classification work. Some 
xvsource categories have very heterogeneous light reflectnnct- patterns 
(liindnat spectral signatures). Some are interspersed among other cate- 
gories. Many wildland vegetation categories arc very heterogeneous e* d 
obviously require special treatment. 


I 
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We will uce "clustering" for this heterogeneous condition. Select 
one training area which representc a land-uoc category of intercut yet 
is very diosccted and mixed up in appearance on the false color cor.jo- 
cite. Bocauce one is leas sure of actual ground location when speci- 
fying thece heterogeneous types, a larger training area is preferred; 
say 200 to 100 pixels. Record the coordinates for this area in the 
some format that the other five areas were recorded in. 

Generate Statistics (5 arena) 

landsat spectral statistics arc the driving force behind the inure 
classification process. They are us * oy the maximum likelihood clas- 
sifier to calculate the probability of a pixel falling into n particular 
class. It is central to classification accuracy that the«.c statistics 
be us homogeneous as possible (i.e., low variance in each channel) an c 
yet represent the land cover tyj* os accurately as possible. 

We nre tuking the supervised approach to rpectraJ statistics develop- 
ment lor the [/ homogeneous training areas we Iwve selected. We aimp?y 
tally all pixels occurring within the training area and compute Die 
Dulti-variate means, variances and covariances for the data set defined 
by the training area boundaries. 

The VIC Ah program "STATS" is used to generate supervised training 
statistics. A complete VICAR run is explained in Appendix 1. 

Cluster (i area) 

We feel tliat our heterogeneous land-use category probably cor tains 
KVi-ral spectral c lueses which will likely color the sane color (i.c., 
place in the same final class) in the final pnoto product of the 


finished classification. Yet we are unsure of where theee spectral 
classes are located (i.e., which pixels are which) and how aany are 
there. Statistical clustering will allow the Identification of 
natural grouping or cluatering in the data eat. to will uat the 
VI CAS prograa "CLUSTER" to autoaat iraily eelec* 10 spectral classes 
Iron our heterogeneous area. Figure l for details oo using 

"CLUSTER." 

Cluster Entire Training Area 

It is often useful to cluster the entire stud/ area in order to 
pick up spectral statistics which as/ not be in the training areas you 
have defined. We will not take the tine to do this additional clus- 


tering during the workshop 
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VI. Merge and Edit Statistics 

An analjat cannot assuae that the statist ica produced by the super- 
vised approach and by clustering are statistically distinct fros each 
other since the work was done independently. A classification of pas- 
ture produced fros clustering nay have near identical spectral patterns 
as a pasture class developed fros a pasture training area. Therefore, 
it is not enough to sisply coabine all of the statistics produced fros 
various approaches; they sust be edited. 

It is quite possible that two classes say be sinilar enough spec- 
trally to cause ais-classification. When two spectral classes overlap 
exceasively in spectral space, the computer experiences difficulty in 
deciding which class to place pixels into. One tool which assesses the 
spectral siailarity or uniqueness of two classes is the separability 
statistic. Separability statistics are computed for each pair of spec- 
tral classes and coabined in a separability aatrix. In this exaaple, 
two classes were very sinilar with a separability of 0.13* These 
classes, 4 and 5. were developed fron the two cypress training areas 
which seeaed to appear slightly different on the red (channel 2) line 
printer aap. This low separability statistic indicates that these two 
classes could not be accurately distinguished by the classification 
algorithm. Since neither of these classes conflicted with other classes, 
they should be pooled together to fora one spectral class for cypress 
swamp. This was done with the following core VICAR coanands: 

E,STATEDIT,IN f GKP, ,PS1 

P,PS1 

SE1JECT , 1 ,2,3*6, V, 8,9, 10, 11 , 12,13* 1^* 15* POOL,4,5 
SEPAR 


:;uO?« Vw MGE 

f ip T*ook Ti 
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A two channel map is also uasful for visually comprehending the 
positions of each class in spectral space. A two channel plot is 
drawn in Figure 1. Classes 1-5 were developed by a supervised approach 
and show wide dispersion except for k and 5* This is to be expected 
since we trained on vastly different resource categories to fora 
classes 1-3* Classes 6-15 were developed froa clustering the hetero- 
geneous mountainous area near the urban center. Note how the clus- 
tering has produced a "cloud” of similar yet spectrally distinct 
classes. This is typical for classes developed froa clustering. Also 
note that the euclidian distance between aany of the unsupervised clus- 
ters is as aaall as the distance between class and 5* The variances 
are smaller for the clustered data, however, allowing for less spec- 
tral overlap. In general, clustering allows the identification of 
more 6pectral ''lasses than the supervised approach. 

We did not experience any spectral overlap between classes developed 
from clustering and supervised classes. We were lucky. In a more 
thorough development of a classification, one would have trained on more 
upland, mountainous areas and experienced some spectral overlap and 
confusion. When the editing process becomes more complex, two good 
general rules to follow are these: 1) When two classes are very similar 

yet don't overlap with other classes, they should be pooled. 2) When a 
class is similar to two or more otter classes the class causing the prob- 
lem should be deleted. Of course many other more subjective considera- 
tions affect the decision to pool or delete, such as the number of pixels 
in the classes, the importance of a particular class to the user and the 
variances of the classed involved. 
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VII. Final Classification 

The VICAR software has two classification algorithms available, 

HAILS and FABTC1AS. RAILS is a suit i spectral classifier which uses 
the Bayesian naxiaua likelihood algoritha. FASTCLAS is a sultispectral 
classifier using an algorithm which cosbines the parallepiped and 
Bayesian techniques. The classification is done in a two-step process. 
First, parallepiped decision boundaries are set up at plus-or-sinus a 
nuaber of standard deviations froa the aean in each channel. The user 
specifies the nuaber of standard deviations to use. In order for a pixel 
to be placed into a specific class, the aean value aust fall within the 
decision boundaries for that class. If a pixel's spectral signature 
falls outside the decision boundaries for all classes, the pixel is 
assigned to the unknown class. A pixel whose spectral signature falls 
within the decision boundary for aore than one class is considered ambi- 
guous. The user can resolve the aabiguity by classifying this pixel wit: 
u maximum likelihood algorithm. The following basic VICAR commands were 
used for FASTCIAS: 

E, FASTCLAS, (A, GRP), CAT, ,PS2 

P.PS*’ 

MSS, 4 SIGMA, ’>.0 

These statements will execute FASTCLAS using an input data set, 

"A," and a statistics file, "GRP." The output file is called "CAT." 

A complete listing of the VICAR command sequence to edit the statis- 
tics and run HAYES and FASTCLAS is presented as Appendix II. 



APPENDIX 1 


Annotated VICAR Cummamd Saquance to; (1) Generate Statistic a. 
(2) Cluster to Create Statiatica and (3) Edit Statiatica. 


Each section of the following VICAR run ia lettered and the explanation for 
each section is located after the letter aa a footnoto. 


Line # 



SETUP, NEED TAPE VJT 008 9 TR800 BPI WITHOUT RING 
TI*,500 

READ,* ,WJT008,N,9 

SAVE , 2 , 101 , 101 , VICAR1 ,COLST .COLCLUSM 
RESERVE , 4 , 514 , 5 lA , • , ( C0L1 , C0L2 , C0L3 , COIA ) 
RESERVE, 1,2055.514,*, A 
E,SAR,*N0l,C0L1,( 21 1,328, 512, 512) 

E ,SAR . *NQ2,C0L? , (21 1 , 328 , 512 , 512) 

E, SAB, *1103, C0I<3, (21 1,328, 512, 5 1 2) 

E , SAB, •H04.C01A, (211,328,512,512) 

E .STATS ,( C0L1 , C0L2 , COW , COIA ), C01ST ,, PS 1 
P,F 61 

HIST, SPEC, TRAIN 
CLASS, 1,240,430,20,20 

t 1 

CLASS.5,466,853,20,20 
E ,MSS ,(C0L1 , C0L2 , C0L3 , COlA ) , A 
E, CLUSTER, A, COLC LUSH, ,PS2 
P,PS2 

MSS ,4 srec NCATS .10 NCMIN , 5 
HIST, 2 , 4 DELTA ,0.45 CONV.tf.O 
SAMPLE, 4 TRAIN 
cuss , , 328 , 500 , 100 , 100 
E .STATED1T ,COIST , • , ,PS3 
E .STATED1T , COLCLUSM , • , ,PS 3 
P.PS3 

COPY .ALL CHECK SEPAR 

END 


EXPLANATION OF LETTERS 


A. The:*.<‘ statements are standard set up cards for a VICAR run. 

Lane 1 , tolls the operator what tape to use 
Line 2, sets an upper time limit of 500 seconds 
Lane 5, reads the tape and labels the tape "N" 

Line 4, reserves 2 permanent data sets called "OOLST" and "COLCLUSM" on a 
disk pack called VICARl. 

101 samples per line by 101 lines 

Line 5» reserves 4 temporary data sets in the public memory (*) called 

C0L1, COL?, C0L3, C0L4. Each data aet is 514 samples x 514 lines. 
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Line 6, rcccrvec one temporary uatt set, 2055 samples by 514 lines in the 
public memory. Tho file if? to be called "A." 

Lines 7-10, request that four hands of Inndsat data be copied fron taj*? "11" 
to disk files called C0L1-4. (COL for Columbia) Only a subwindow 
of data is to be copied however, storting nt line 211, sat pie .32° 
and continuing 512 samples by 512 lines. 

B. These ctutcment6 generate spectral statistics for 5 training areas. 

Lines 11-3 0, generate spectral statistics using 4 hands of londsut data 

(C0L1-4) and call the output file C0L5T (for Columbia statistics). 
The sine field is defaulted (,,) and the statistics are to be 
generated according to parameter statement 1. Parameter statement 
1 (lines 12- IS) calls for histograms of all input bands of each 
spectral class and spectral plots of all bandb. The word "TRAIN" 
tells the computer that training fields are to follow. Lines 14-18 
give rectangular window coordinates (starting line, starting sample, 
number of lines, number of samples) for 5 training areas. 

C. This one statement (line 19) reformats the 4 bonds of data into an KSS 

format window called "A." 

D. These statements request a clustering algorithm to be passed through data 

set "A." 

Line 20 , requests the "CLUSTilR" routine to he used on data set A. The 

output file is to be called COLCLUCK (for Columbia clustering on 
a mountainous area) and the cize field is defaulted (,,). The 
parameter statements (lines 22-25) indicate that; (l) the input 
data is 4 channel, MSS format, (2) spectral plots of all channels 
arc* requested, (3) number of classes to bo generated is 10, (4) 
minimum number of classes after merging 1g 5. (5) histograms are 
requested of channels 2 and 4, (6) the delta or minimum separability 
statistic is 0.45, (7) the poicent convergence required is 
(8) every 4th pixel is to be sampled and (9) one training area is 
to be submitted, 100 lines by 100 samples. 

K. This section requests that separability matrices be printed out for both 
statistics filo6 created, C0I£T onu COLCUJSM. 


t 
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APPENDIX 2 

VI CAM coMand sequence to adit the atatieticn end run BAYES 
and FASTCLAS. 


CHARGE ,*4352 
TIME, 1900 

SETUP, NEED TAPE WJT0f)8 9TR ttOO BPI WITHOUT KING 
SETUP, NEED TAPE WJT018 9TB 8(X) BPI WITH RING 
HEAD,",WJT0O8,N,9 
WBITl,»,WJTOl8,T,9 

RESERVE , 4 , i>l4 , 514 , • , ( COLl , 0012 , 0013 , COLA ) 

RESERVE , 1 , 2055 , 51** , * , IN 
FIND , (QRP, VICAJM ) 

E, SAB, *1*01,001/1, (211, 328, 512, 512) 

E,SAR,*NQ2, OOL2, (211, 328, 512. 512) 

E,SAR,*N03, 0013,(211 ,328,512,512) 

E ,SAB , *N0<» , OOIA , ( 211 , 328 , 512 , 512 ) 

E , NSS , ( OOL1 , 0012 , 00D3 , OOIA ) , IN 
E,FASTCIAS,(IN ,GRP),00L1, ,PS2 
P.PS2 

MSS, 4 SIGMA, 9.0 

E .STRETCH ,C0L1 , C012 ,,( LINEAR ,0 , 63 ) 

E, DISPLAY, 0012,*, ,PS3 
P.PS3 

• ,-=+/W*0W=)(4S ' 

• +N/ * / 

• • • 


5,5 

E .BAYES , ( IN ,QRP) ,0013 , , (MSS ,4 ) 

E, STRETCH ,0013,00 LA, .(LINEAR, 0,63) 

E, DISPLAY, 00IA,*,,PS4 
P,PS4 

• .-=*/V*0W=)(4S 

' *H/ = / 

• • • 

5,5 

E,SAR,OOL1,T 
E, SAB, 0013, T 
END 
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Project Summary 


In general, the module was well received by the VICAR users and 
potential VICAR users at the NASA Workshop. There was enough 
time allotted for the lecture and exercises. The material seemed 
to be addressing the appropriate level of sophistication for 
efficient technology transfer. 

The schedule called for the students to actually run VICAR soft- 
ware, "live" at the workshop. This was quite beneficial for 
establishing user confidence, yet required a tremendous amount of 
backup effort. Mike MacDonald of Technicolor gave his energy 
above and beyond the call of duty. If Mike had not been willing 
to put in three 18 hour days, the "live" portion would have 
failed miserably. 

The lesson learned was not to underestteate the amount of support 
required for students In a workshop to actually process Landsat 
data. I would suggest that 4 or 5 consultants be available to 
review student card decks before submission to the computer and 
to help de-bug VICAR runs after they have failed to operate. 

The discussion session would have been improved substantially If 
more examples of VICAR classifications hud been available. The 
Principal Investigator came up short in this area. In the future, 
it would be best to have several comparative VICAR classifications 
(o.g., KASTCLAS and BAYES with supervised statistics, with guided 
clustering and witli unsupervised statistics) available for student 
review and discussion. 1 was counting on the student output which 
didn't materialize to the level required for an in-depth discussion. 


SECTION TWO 


FOREST RESOURCE ASSESSMENT 
OF 

HUMBOLDT COUNTY 
By: 

Lawrence Fox III and Kenneth E. Mayer 
SCHOOL OF NATURAL RESOURCES, DEPARTMENT OF FORESTRY 
HUNSOLVT STATE UNIVERSITY 


A Report of Work Performed Under 
NASA Grant No. NSG-2341 


NASA ’ '“l CAL MONITOR: 

Su.. M. Norman 


NASA TECHNICAL CONSULTANTS: 
James Brass and David Peterson 


In Cooperstion With 

FOKfcSr RL SOURCE ASSESSMENT PROGRAM (FRAP) 
and 

CALIFORNIA STATE DEPARTMENT OF FORESTRY (CDF) 


FRAP COORDINATOR: 
Nancy Tosta-Miller 


INTRODUCTION 


On January b, i960. Judge Lawrence K. Karlton of the U.5. District 
Court of Eastern California ruled that the National Environmental Impact 
Statement (NEIS) supporting the second Roadless Area Review (RARE II) did 
not satisfy the requirements of the National Environmental Policy Act 
(NEPA). The court then enjoined the U.S. Forest Service (USFS) from 
developing any of the 47 disputed non-wilderness areas in California 
pending proper consideration of wilderness values In compliance with 
NLPA. Primary deficiencies cited by the judye Included the "lack of 
site-specific data on RARE II areas" (J. of Forestry, 1980). 

This dramatic court decision with its far-reaching national implica- 
tions hi yh lights a pressing need for Inventory data on a site specific 
basis. Landsat, a NASA-sponsored earth imaging satellite. Is providing 
the potential for site specific forest inventories of very large areas. 

One Landsat scene covers approximately 3.4 million ha (8.46 million ac) 
or 34,255 km? (13,214 square miles) with t. grid cell Information base, 
each cell having a ground area of approximately 0.46 ha (1.12 ac). The 
satellite measures reflected "light" In four bands of the electromagnetic 
spectrum, two visible and two non-vislble Infrared wavelength bands. 
Reflectance patterns called "spectral signatures" are measured by the 
satellite as it passes in sun synchronous, near polar orbit. The space- 
craft returns to the same location on the earth every 18 days (Sabins, 
19 / 8 ). 

The legislators of California have responded to forest inventory 
needs by passiny the California Forest Resources Assessment and Policy 
Act of 19/7 (FRAPA). This bill (AB 452) established a project to improve 
the information base, upon which forest policy decisions are formulated. 
Under policy guidance from the California Board of Forestry and with the 
help of the Secretary of Resources, the Director of Forestry was required 
to prepare a forest resources assessment and analysis by July 1, 1979. 

The bill also provides for an updated assessment every five years. 

The staff of tne Forest Resource Assessment Program (FRAP), a branch 
of the California Department of Fuiestry (CDF), cooperated with the NASA 
Ames Research Center (NASA-ARC) Moffett Field, California, to use Landsat 
data in the preliminary assessments. Forest cover was classified by 
computer, using Landsat data of the entire state. Land condition and type 
of vegetative cover was placed into one of 17 categories which included 
conifer forest, coni fer/hardwooa forest, hardwood forest and brush. 

Area summaries were provided by category for each county of the state 
(Peterson ana Tosta, 19/9). 

To further develop the detail of tne assessment program, a cooperative 
project was undertaken in 19/9-80 between CDF, NASA, Humboldt State 
university (HSU), California Polytechnic State University at San Luis 
Ubispo, and the Tahoe National Forest. The objective was to intensively 
inventory three areas in California with Landsat digital data to determine 
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it more classification detail could be achieved on a local basis than 
had been possible in the state-wide survey. The results would assist 
CDF in fulfil liny future forest assessment responsibilities. 

This report is a description of the Landsat Inventory of Hunboldt 
County, completed by HSU. It therefore represents one phase of the 
cooperative project. The goals of the Humboldt project were to: 

l) Develop an intensive forest land Inventory strategy for 
Humboldt County that could evaluate the full discriminatory 
power of Landsat spectral signatures for vegetation classi- 
fication. 

i) Map and inventory the forest land of Humboldt County using 
Landsat digital data via computer classification techniques. 
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PROJECT AREA 

Huaboldt Is California's 14th largest county comprising 9254 km 2 
(3573 square miles) (Figure 1). The county Is approximately 170 km 
(105 miles) long, north and south, and 56 km (35 miles) wide, east 
and west. Humboldt Is primarily mountainous with the exception of the 
alluvial lowlands In the Eel River Valley and Humboldt Bay Area. Eleva- 
tion ranges from sea level to over 2133 m (7000 feet). Precipitation 
occurs during the winter months, with an average rainfall of 101 cm 
(40 Inches) with excesses of 229 cm (90 Inches) occurring In Isolated 
coastal mountain areas. Temperatures are highly variable In the county 
with lows for the winter of -8.8°C. (16°F.) In the mountainous areas 
and 12.2°C. (54°F.) In the Eureka area. Sumner high temperatures range 
from 18°C. (60°F.) on the coast to 40°C. (100°F.) In the Inland areas. 

The 1980 census showec 107,500 county residents with approximately 
50 percent living within a :5-m11e radius of Eureka. Humboldt has a 
higher population than any of its neighboring counties, Del Norte to the 
north, Mendocino to the south, and Trinity and Siskiyou to the east. 

Twenty-seven percent of the county Is In public (federal, state, 
county, and municipal) ownership and 73 percent Is private (Himiboldt 
Atlas, 1975). The Forest Products Industry provides approximately one 
job in four. Humboldt is the top lumber- producing county In California 
and ranks among the top lumber-producing counties In the United States. 
Redwood (Se quoia se mpervi rens ) Is the major commercial species with 
Uouglas-fl r ( ftseudotsuqa menzlesi 1 ) and hardwoods providing the remainder 
of the important timber resources. Old growth timber Is limited, creating 
a major shift towards second growth harvesting and management. 

Important conifer species are: redwood, Douglas-flr, grand fir 

(Abies gr andls ) , sitka spruce ( Plce a sltchensis ) , Western hemlock ( Tsuqa 
hete ropnylla ) , Port Orford-cedar ( CTTamaecyparls lawsoniana ) and Western 
re'Jcedar ( Thuj a plicata ). Important hardwood species are: tanoak 

( Li thocarpus densiporus ) , Pacific madrone ( Arbutus menziesil) , giant 
chinquapin ( Chryso T epis c hrvsophyl la ) , red alder ( Alnus rubra ) , Oregon 
white oak (Quercus garryana ), and California black oak ( Quercus kellogqii ). 
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INVENTORY DESIGN 

Landsat provided a comprehensive data base with a grid cell 
coordinate system. The Inventory was designed to classify every picture 
element (grid cell) and therefore was a 100 percent Inventory with no 
sampling. While It was certain that all of the area would be Inventoried 
to a .4? ha cell size, the exact classification was not prepared before- 
hand. The spectral signatures developed from the Landsat pixels may nave 
represented various levels of classification detail depending on environ- 
mental conditions, forest types Involved, the discriminatory power of the 
computer programs, and the skill and experience of the analyst. 

We began the project with the goal of Identifying two major conifer 
species types, redwood and Douglas-flr, and various mixes, as well as 
several forest density categories and size classes. A literature review 
of work completed prior to 1975 Indicated that we would not reach this 
level of detail (deStinger, 1978) . Our experience had been that the 
Initial classification goals were achievable (Fox and Mayer, 1980) and 
recent work by others Indicated a high probability of success (Walsh, 1980). 


ANALYSIS TECHNIQUES 

The classification of Lanosat data can be divided Into three dlctlnct 
tasks: 

1) Developing a set of descriptive statistics for the spectral 
classes required. 

2) Classifying the study area Into these classes. 

3) Identifying the resource label for those spectral classes. 

To complete this project, we used unsupervised and guided clustering 
techniques (Flemming, 1975) to develop the needed spectral classes and 
multivariate discriminant analysis for classifying every pixel In the 
study area (Sabins, 1978). 


Scene Selection and Image Mosaic 

After a review of available Landsat Images of Huntooldt County that 
were cloud- free and date sequential, two April 12, 1977 scenes were 
selected. Each image covered approximately SO percent of the county. 

A computer compatible tape was acquired from Earth Resources Observation 
Systems (EROS) Data Center. Radiometric anomalies including bad data 
lines and points, radiometric striping and atmospheric scattering were 
not corrected or normalized before processing. These errors were not 
significant for the data tape. 


The two Images were overlapping north to south by approximately 
200 lines. To assure an accurate line-by-line natch for joining the 
two Images, the overlapping lines were removed, which was completed on 
the IBM- 360 computer at NASA-ARC. This assured continuous, non-overlapping 
data north to south. Each scene was copied to tape as well as to computer 
disk for easy access during processing. 


Resource Category Selection 

In order to determine the maximum detail that could be obtained from 
Landsat digital data classifications, resouce categories were selected 
for targeted forest land variables. The approach focused on conifer species 
differentiation, tree size class, and crown closure category* Furthermore, 
forest categories such as hardwood, brush, and regeneration areas were also 
being selected for classification. 

Before a feasible classification scheme was selected, an examination 
of tne vegetation communities was completed. As the different forest 
comnunltles were Identified, a determination was made relative to the 
following criteria: 

1) Is the resource type an Important part of the vegetation 
mosaic of Humboldt County? 

2) Is there enough acreage of that vegetative community to 
facilitate computer processing 7 

3) Would the resource type be spectrally separable from the 
other types? 

If It was determined that a vegetative community met these criteria. It was 
Included Into the preliminary classification scheme. The CAlVEG classifi- 
cation scheme was used as reference to develop the vegetation types (Parker 
and Matyas, 1980). It was also decided that resource categories and 
classification goals would change throughout the class! ficatlon process as 
the limitations of spectral discrimination and computer software were 
encountered. 


Training Area Selection and Groundtruthlng Effort 

Approximately 10 homogenous training fields were selected for each 
resource category. The training areas were at least 10 ha (25 ac) In 
size, with the majority of training fields being »20 ha (50 ac). The 
training areas were located on U-2 1:32,300 color Infrared aerial photo- 
graphy and outlined on U.S. Geological Survey (USGS) 7V quadrangle maps 
to facilitate groundtruthlng. 

A field visit was made to each training area to determine plant 
species, vegetative crown closure, age, size class, aspect and slope. 
TMs information was recorded and used as an aid In the labeling of 
spectral classes generatcJ. 


Landsat Classification Procedures 


In order to extract the "raw" (unaltered) Landsat spectrel Informa- 
tion from within each training field, the Landsat Imaaes were displayed 
on the Interactive Digital Image Manipulation System (IDIMS) cathode ray 
tube. The training areas were located by visual Inspection of the Landsat 
data. These areas were removed from the data set and copied on disk to 
facilitate guided clustering. 


Guided clustering was completed or. each set of training areas. EDITOR 
software was used to analyze the Landsat data. Procedures used were similar 
to those reported by Fox and Mayer, 1979; Forbes et al. 1980; and Mayer et 
al. 1980. 


Unsupervised clustering was performed on each Landsat Image to 
develop spectral statistics unique to the resources present in the scenes. 
The two statistics files were compared for spectral separability (Fox and 
Mayer, 1979). The statistics files were edited to eliminate classes with 
spectral similarity. The final statistics were merged with the statistics 
developed from guided clustering. 


The final, edited statistics file was used as input to a maximum 
likelihood decision function. This step completed a multivariate 
discriminant analysis of the Landsat data and classified every pixel 
into one of the spectral classes developed. 


Many more spectral classes (unique light reflectance pattern) are 
usually generated than tne number of resource categories defined (Mayer 
and Fox, 1979) . in order to comprehend the data, several spectral classes 
must be grouped together into resource categories. Spectral classes were 
grouped according to geographical proximity and similar resource classifi- 
cation as observed on color infrared photography. Spectral curves of 
Landsat generated spectral classes were also used to determine the proper 
resource identity (Mayer and Fox, 1980). 
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Accuracy Evaluation 

The Landsat classified scenes (April 12, 1977) were evaluated using 
U-2 color Infrared photography (Nay 8, 1978). Classified data windows 
were selected from both scenes that could be easily located on the ground 
and on the photographs. Single pixels were randomly selected from these 
windows. A two man team was used to complete the evaluation. One member 
of the team located the pixel to be evaluated on the photographs. The 
second member completed the photointerpretation to determine the "truth" 
(actual resource Identity). 


by: 


The overall mean probability of correct classification was calculated 


P 


Y 

7 


LJHEKt: Y Is equal to the number of correct Landsat pixels 

and X Is the total number of pixels Identified by 
photol nterpretatl on . 

For a complete technical description of the evaluation process, refer to 
Mayer and Fox, 1980 and Fox et al. 1981. 


Area Summary 

Area summaries for the Landsat resource categories were produced 
for Humboldt County. The boundary of the county was digitized Into the 
computer data base in order to create a specific geographical segment 
of the categorized data. Area summaries within that segment were accessed 
through the EDITOR software system. 
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RESULTS 

Resource Category Description 


The Landsat spectral data was classified Into 64 classes, each one 
representing a statistically unique light reflectance pattern (Table 1). 
The spectral classes were assigned to 15 resource categories based on 
similarities In composition (Table 2). The forested land was classified 
into nine categories, based on species composition (redwood through hard- 
wood). All of the forest categories represented closed canopy forest 
stands with various proportions of hardwoods and conifers. The total 
vegetative crown closure was consistently greater than 65% for all nine 
forest categories. 

The "Redwood" category Is dominated by redwood ( Sequoia sempervirens) 
which represents 80% of the conifer trees present. Up to 20% of the 
conifers present may be one or more of the following species: Douglas-flr 
( Pseudotsug a inenzlesii ). grand fir ( Abies grandls) , sltka spruce ( Plcea 
si tchensisT T and western hemlock ( Tsuga heterophylla) . 

The "Redwood" class may contain small amounts of hardwood, usually 
red alder ( Alnus rubra ). The stands are deflnetly conifer dominated; 
however, with hardwood seldom occupying more than ten percent of the 
forest canopy. 

The "Uouglas-fl r ^>80%" category also represents a closed canopy 
conifer stand. Douglas-fir dominates and comprises over 80% of the total 
conifer cover. Associated conifers may occupy up to 20 % of the stand and 
be any or all of the following species: ponderosa pine ( Pinus ponderosa ) 

sugar pine ( Pinus lambertianna ) , incense-cedar ( Libocedrus decurrens ) , 
white fir ( Abies c oncolor ), c?d noble fir ( Abies~procera) 7 Small amounts 
of hardwood (less than lOx of total trees) may also be present In the 
stand. 

The "Douglas- fir^ardwood" categories (61-80% and 26-60%) are fully 
stocked stands containing a mixture of Douglas-flr and hardwood. Douglas- 
fir dominates, representing 61 to 80% of the trees present In the 61-80% 
category. Hardwoods are principally tanoak ( Lithocarpus denslflorus ) and 
Pacific madrone ( Arbutus menziesn ) . Hardwoods dominate the 2(J-6&% 
category. Since Douglas-flr accounts for 26-60% of the trees present, 
hardwoods and/or brush comprise the remaining proportion of the stand. 

Up to 20% of the conifers In these stands may be comprised of species 
other than Douglas-fir, the same associates listed for the previous 
Douglas-fir class. 
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Average stand height of dominants and co-dominants. 

Total Hardwood crown closure. Other vegetation always present to 
make total vegetative crown closures consistently above 65 %. 
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The "Mixed Conifer" category represents a closed crown stand In 
which no one conifer species dominates. The class Is made up of recfcood 
and Douglas-flr near the coast and pine, fir, Dougl«s-f1r, and occasionally 
white fir further Inland. 

The "Mixed Conifer/Hardwood" classes are Identical to the Douglas-flr/ 
Hardwood classes discussed previously except that the conifers In these 
categories are mixed with no single species previously predominating. The 
coastal forest community described as "redwood" was placed In this category 
(Parker and Matyas, 1980). These forests often contain 30 to 50% hardwood 
and the conifers rarely are redwood only. Often 30% of the conifers present 
are other than redwood. The species listed In the "redwood" category are 
the major ones present. 

The two "Hardwood" classes are comprised of red alder, tanoak and 
madrone In various densities. These classes contain woody vegetation 
greater than 15 feet In height. This Is an Important distinction since 
tanoak grows In brush form as well as tree form. Often the only difference 
between brush and hardwood Is a difference In canopy height. The Harfrood, 
>60% crown closure class contains greater than 60% hardwood forest canopy 
closure. The understory Is brush or grass/forest. The HarAvood, 26-60% 
crown closure class contains between 26 and 60% crown closure of hankvoods. 
Brush or grass/forest make up the remainder of the total vegetative cover. 

The "Soil" category Is primarily bare soil yet may contain up to 
25- vegetation, with scattered trees, brush or grass/forest. 

"Agriculture/Grasslands" is a combination range/agriculture category 
which contains both Irrigated and non-lrrlgated pasture. 

"Brush" is a category describing woody vegetation less than 15 feet in 
heiyht. Tanoak, blue blossom ( Ceanothus thyrsiflorus ) , snow brush ( Ceanothus 
velutinus ) , deer brush ( Ceanothus integerrimu s), and coyote brush ( Baccharis 
pilutaris) are major species (refer to Parker and Matyas, 1980 for a 
complete listing). The brush category may contain conifer seedlings that 
are beneath the brush canopy level. This class does not contain significant 
amounts of conifer seedlings or saplinys protruding above the brush canopy. 

The "Brush/Regeneration" category Is identical to the brush canopy except 
that it contains significant amounts of conifer seedlings and/or saplings 
protruding above the brush canopy. The conifers represent over 40% of the 
canopy in this class. Conifer regeneration may be any of the conifer species 
described previously. 

The "Other" category describes lands of little importance to forest 
management, urban, gravel and snow. 

The "Water" category describes rivers and marine estuaries in the 
county. 
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Accuracy Evaluation 

The overall mean probability for correct classification was 8?. 2 
percent. The accuracy estimate was determined from evaluating 320 single 
pixels scattered randomly throughout the county. Individual category 
accuracy statements were not generated as photo coverage was not complete, 
creating problems in finding and locating adequate sampling areas. Further 
evaluation of detailed resource categories will take place In the future 
(Prime Timberland Contract 3-080-078), as photographs become available. 


Area Sunvna ry 

Landsat identified 962,291 ha (2,309,512 ac) In Hunboldt County 
(Table 3). It was found that Douglas-fir was the most predominant conifer 
category, comprising 21 percent (485,481 ac) of the total land base. This 
was expected, as a preliminary survey of the county suggested that there 
was more Douglas-fir present than redwood (Plank, 1974 ). 

Reckvood>80% comprises four percent of the total land base; approxi- 
mately 37,928 ha. It is important to note that there is more redwood in 
the county than is revealed by this figure. Many of the "second growth" 
redwood stands were placed into the mixed conifer category; primarily due 
to the heavy intermingling (>>20“.) of associated conifer species. In the 
past, many timber inventories referred to any stand containing redwood as 
a redwood stand. Since Landsat spectral signatures respond to total 
vegetative cover, it was not practical or feasible to misinterpret these 
stands as redwood. Consequently, mixed conifer, >80% cc, was the second 
most abundant category, comprising 13 percent (125,300 ha) of the land 
base. 


Conifer forest as a whole makes up more than 60 percent (580,669 ha) 
of the land base 'n the county. This is a graphic example of the wealth 
of timber resources in Humboldt. Furthermore, pure hardwood comprises just 
over 6 percent (62,934 ha). This total may seem to be low, but actually 
it is very realistic. Much of the hardwood in the county does not exist 
in pure stands, devoid of conifers. Hardwoods are found in much of the 
Uouglas-fi r/Hardwood 61-80 percent, 26-61 percent, Mixed Conifer/Hardwood 
ol -80 percent, and 2b-60 percent. These categories comprise approximately 
22 percent (215,158 ha) of the total land base, a significant proportion 
of the forest resources in the county. 

Landsat identified 89,947 ha of brush and 96,198 ha of brush- 
regeneration , approximately nine and 10 percent respectively. Over 
half of the brush lands identified were indicated as being overtopped 
with conifers. Landsat data analysis at this level provided no indica- 
tion of the age structure of the regeneration zones. 


Table 3. Area Estimates For Humboldt County, California Derived 
From Classified Landsat Data - April 12, 1977. 
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RESOURCE LABEL 

ACRES 

HECTARES 

PERCENT 

Redwood >80XCC* 

91,029 

37,928 

4.0 

Douglas-flr >80SCC 

485,481 

202,283 

21.0 

Dougl as- f 1 r/ilardwood 
61 -80'tCC 

55,116 

22,965 

2.4 

Dougl as- fir/Hardwood 
26 -61 ICC 

50,629 

21,095 

2.2 

Mixed Conifer >80%CC 

300,722 

125,300 

13.0 

Mixed Conifer/Hardwood 
61-80XCC 

214,109 

89,212 

9.3 

Mixed Conifer/Hardwood 
26-61 CC 

196.528 

61,886 

8.5 

Hardwood >60 vCC 

74,648 

31,103 

3.2 

hardwood 26-60.CC 

76,396 

31,831 

3.3 

brush 

215,873 

89,947 

9.4 

Brush -Regene rat ion 

230,876 

96.198 

10.0 

Agri cul ture-Grasslar,ds 

152.608 

63,587 

6.6 

Soil< 20XCC-VEG 

93,451 

38,938 

4.0 

Other 

55,146 

22,977 

2.4 

Water 

16,900 

7,041 

0.7 

TOTAL 

- 2, 109,512 --- 

— 962,291 — 

— 100.0 
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Approximately 63,587 ha (152,608 ac) of the Agriculture/Grasslands 
category were found in the county. This total is somewhat lower than the 
information reported from the Landsat statewide Inventory (approximately 
400,000 ac). Much of the rangeland included as grassland previously Is 
possibly being classed as brush In this effort. 
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CONCLUDING REMARKS 


The project goals have been met, as evidenced by the documentation 
included here and the fact that CDF has continued the work begun 
here through an Inter-Agency Agreement with Humboldt State Univeraity. 

The CDF is moving forward to incorporate Landsat data Into their 
operational forest Inventory procedure. 

It is unfortunate that the U.S. Forest Service chose not to Incorporate 
Landsat at this time. We continue to work with the McCloud Ranger 
District on an Informal basis and look forward to future Landsat projects 

Please receive my hearty thank you for the opportunity to work in this 
exciting aspect of modem day remote sensing. 



